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Abstract 

We perform a general analytic study of feasibility of obtaining a com- 
bined explanation for the deficits in the solar and the atmospheric neutrino 
fluxes with two large mixing angles in supersymmetric model with bilinear 
R parity violations. The required hierarchy among the solar and atmo- 
spheric neutrino mass scales follows in this framework in the presence of 
an approximate Higgs - slepton universality at the weak scale. The solar 
mixing angle is shown to be related to non-universality in slepton mass 
terms specifically to differences in soft parameters of the first two leptonic 
generations. It is shown that this flavour universality violation should be 
as strong as the Higgs-slepton universality violation if solar neutrino mix- 
ing angle is to be large. The standard supergravity models with universal 
boundary conditions at a high scale lead to the required Higgs-slepton 
universality violations but the predicted violation of flavour universality 
among the first two generations is much smaller than required. This model 
therefore cannot provide an explanation of large solar neutrino mixing an- 
gle unless some universality violations in soft supersymmetry breaking 
parameters are introduced at a high scale itself. 

1 Introduction 

Experimental observations of deficits in the solar and atmospheric || neu- 
trino fluxes have provided concrete ground to believe in neutrino oscillations. 
These experimental results are consistent with a simple picture of three active 
neutrinos mixing with each other. Within this picture, two independent (mass) 2 
'presently at INFN, Sezione di Padova, via F. Marzolo 8, 35131 Padova, Italy. 
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differences (A , A atm ) among three neutrinos govern the oscillations of the so- 
lar and atmospheric neutrinos respectively. One needs A / A atm < 10~ 2 . Two 
of the mixing angles determining amplitudes of these oscillations are required 
to be large ||. The third mixing angle measured by the survival probability of 
the electron neutrinos in laboratory experiments such as CHOOZ is found to 
be much smaller < 0.1 

Different theoretical possibilities have been suggested for obtaining the above 
neutrino spectrum with two large mixing angles ||. One potentially interesting 
possibility in this regard is super symmetric standard model containing bilinear 
R parity and lepton number violation |g, 0, ||, ||, I(| 11 , 12, l||, |l4| . The following 
features of the model make it an ideal candidate for the description of neutrino 
masses. (1) The lepton number violations and hence neutrino masses and mix- 
ing are described in this model in terms of three parameters. Ratios of these 
parameters control neutrino mixing which can be naturally large. (2) The mech- 
anism for suppression of neutrino masses compared to other fermion masses is 
automatically built-in for two of the most popular supersymmetry breaking 
scenario namely the minimal supergravity model (mSUGRA) and models with 
gauge mediated supersymmetry breaking (GMSB). Extensive studies of these 
models have been carried out in the literature |7], |8], ||, 10]. Our aim in this 
paper is to discuss under what conditions the bilinear model can lead to two 
large mixing angle among neutrinos. We discuss this issue analytically and in 
the process show that the two scenarios mentioned above cannot lead to two 
large mixing angles although small angle mixing solution to the solar neutrino 
problem is possible^. 

The suppression in neutrino masses in mSUGRA and GMSB arises due to 
equality at a high scale(= Mx) of soft parameters of one of the Higgs fields 
(= H[) with the corresponding parameters of the leptonic doublets having the 
same quantum numbers as H[. Small differences arise in these equal parameters 
at the weak scale due to RG scaling. For example, one finds in case of mSUGRA 

Am 2 = (mj, - m 2 H r) « ^| In ^m 2 susy « 2 • 10~ 3 m 2 susy (1) 



feasibility of only small mixing angle solution was pointed out also in |14| . Our analysis 
considerably differs from theirs. 
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where m?,(i = 1, 2, 3),m 2 H i respectively denote the weak scale values of the soft 

i 1 

SUSY breaking masses of the sneutrino and H[ respectively and msuSY is the 
typical SUSY breaking scale ~ 0(100 GeV) . The hb in the above equation refers 
to the 6-quark Yukawa coupling. The neutrino masses in this model involve 
the above and similar differences among B parameters. The suppression in 
these differences leads to suppression in neutrino masses. Thus the smallness of 
neutrino masses is linked to near universality of the Higgs (H[) and sneutrino 
soft parameters. As we will discuss in this paper, the solar neutrino mixing 
angle is directly linked to flavour universality violation, i.e., to differences in 
sneutrino mass parameters themselves. More specifically, the solar neutrino 
mixing angle involves the parameter 

2 2 

6 = ir\ — A» ( 2 ) 

which is required to be 0(1) implying that the weak scale universality viola- 
tion among different flavours are required to be as strong as the corresponding 
Higgs-slepton universality violations. This is in sharp contrast with the expec- 
tations based on mSUGRA and GMSB where the former violations are mainly 
controlled by the muon Yukawa coupling while the latter by the b or r Yukawa 
couplings. Thus 5 in eq.(||) is of 0(1O -4 ) instead of being one. 

Link between universality violation and large mixing was brought out in the 
numerical study of JO]]. In contrast to their work, our analytical study allows 
us to determine specific pattern of universality violation and also allows us to 
quantify the amount of violation needed to obtain the LMA solution for the 
solar neutrino problem. 

We present our results in the following manner. The next section outlines 
general formalism we adopt and our assumptions. It also contains analytic 
discussion of neutrino mixing and masses in this scheme. The close link between 
large angle solar neutrino solution and flavour violation is emphasised in section 
(3) which also contains results based on numerical analysis. The last section 
contains a summary. Some of the technical aspects relevant to discussions in 
the text are elaborated in the appendices. 
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2 Sources of neutrino masses 



In this section, we derive analytical conditions on the low energy susy param- 
eters in order to have a phenomenologically consistent neutrino spectrum. To 
this extent we do not assume any specific structure of the soft masses. Fur- 
ther assumptions regarding various contributions to neutrino mass spectrum 
within mSUGRA inspired scenarios are discussed as and where required. The 
superpotential takes the following form in our case: 

W = h? j Q i u c j H 2 + hfQid^H[ + + ^H[H 2 + e^i? 2 . (3) 

Without loss of generality, we have chosen above the basis in which the down 
quarks and charged lepton masses are diagonal. The e, t characterise lepton 
number violation in this basis. 

We have the following soft supersymmetry breaking terms at the weak scale: 

Vsoft = rn 2 H[0 \H'®\ 2 + m 2 H o\H 2 1 \ 2 + m%,\v % '\ 2 

- (B^'H'^Hl + c.c) - Bi€i (Vtf 2 ° + c.c.) + . . . (4) 

Note that the above equation refers to soft terms at the weak scale. For sim- 
plicity we have displayed only the terms involving neutral fields in the above 
equation. The following comments are needed in connection with eq.(Q): 
(i) Although we have allowed for arbitrary diagonal sneutrino masses, we have 
not included off-diagonal sneutrino masses in this primed basis since such off- 
diagonal masses are severely constrained by flavour violating processes, e.g. 



fj, e7 |5]. 



(it) V so ft does not contain sneutrino-Higgs mixing terms of the form m~, Vi*H'^ 
although they are allowed by the gauge symmetry. Such terms are not present 
in the minimal supergravity theory at high scale. The renormalization group 
(RG) equations for ^n 2 -,^ given in the appendix, eq.(^Tj) show that these terms 
cannot get generated even at the weak scale if they are not present at high 
scale. Thus it is meaningful to omit these terms. We should emphasise that 
this statement is very specific to the particular basis in which bilinear terms 
are not rotated away from the superpotential until the weak scale and neglect 
of such terms would not be justified in any other basis. In our case, the Vi*H\ 
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term would make its appearance when we go to the basis with no bilinear R 
violating terms in the superpotential at the weak scale. 

The neutrino masses arise from several sources in this model. Discussion of 
these sources becomes transparent if we re-express eq.(^) in the new basis in 
which bilinear terms are rotated away from Wft. 



Hl = ? ' 

U = ^7^ . (5) 
A* 

2 

This basis are simple but are orthonormal only up to O(^). This approxima- 
tion is sufficient for most of our discussions since q are required to be much 
smaller than the typical SUSY scale // in order to reproduce the scale of neu- 
trino masses correctly. Generalisation of eq.(||) valid to higher order in e« and 
its consequences are discussed in the appendix B. Eq.(||) takes the following 
form in the unprimed basis: 

W = h^QiU^H 2 + hiQid^Hx + hlL^Hx - X'^LiQjdt - X^UL^l + nH x H 2 . 

(6) 

where 

2 12 , 2 , 2 , 2 ^ 12 

H = [i + e 1 + e 2 + e 3 « fi , 

A ijk / n j jk j 

Xijk = (Sikhf^-Sj^) . (7) 

Similarly, after rotating primed terms in eq.(Q) and adding the contribution 
of the supersymmetric part, we get the following expression for the full scalar 
potential in the unprimed basis: 

V sca l a r = (m^o + ^ 2 )|#iT + (m^O + ^)I# 2 T + + H l + CX 

- (B^HlHl + c.c) - AB^ (yiH% + c.c) 

+ W + 9' 2 )(\h i\ 2 + \^\ 2 -\h 2 \ 2 Y , (8) 



2 Note that this definition of a new basis is same as that of Ref.|n|. However in the present 
work, this rotation is done only at the weak scale in contrast to where it is scale-dependent. 
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where 

Am 2 = m\, - m 2 H ,« AB t = B, - B^ . (9) 

Two major sources of neutrino masses arise from eqs.(|7|,^). Minimization of 
eq.(|8|) generates sneutrino vev: 

< Vi >= tih (10) 

where 

vi (-Am] + tan/3 /xABj) 
l ~ n (m| + l/2M§cos 2(5) [ ' 

v\ =< H® > and Mz represents the Z boson mass. Sneutrino vevs lead to 

neutrino masses through their mixing with neutral-gauginos: 

Mtree = A < Vi X Vj >= A^jkikj . (12) 

Aq is obtained by diagonalizing the 7x7 neutrino-neutralino mass matrix in 
the standard way [17]: 

A = v(g 2 + g' 2 ) (13) 

2 (-cfiM 2 + M$ v sm2f3(c + tan 2 6 W )) ' 1 ' 

where 9\y represents the Weinberg angle and M\y represents the W-boson mass, 
c is given by 5g 2 /3g 2 ~ 0.5 with M2 representing the standard gaugino mass 
parameter. 

The trilinear terms in eq.(|7|) lead to the second contribution to neutrino 
masses at 1-loop level. Since these couplings are proportional to the Yukawa 
couplings, the dominant contributions arise due to exchanges of the 6-quark- 
squark and r-lepton-slepton in the loops. The loop induced mass matrix is of 
the form : 

(Mioop^j = e t ej (A b + A T (1 — 5 l3 ) (1 - 6 j3 )) (14) 



where 



M = hf sin (p b cos (p b ln [t^T ) > ( 15 ) 



1 y i ^3 • j. j. 7 / M 2 
— _ M m</> r co S ^l^ 



A t = TEZaT^r sin0 r coscj) T ln ( ] . (16) 
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Here (j>b,M denotes mixing between the left and the right handed squark (sneu- 
trino) fields. These mixing angles are proportional to the b and t Yukawa 
couplings. Approximating them by ™ b ' T we get the following numerical values 

A Q « 5-10- 3 GeV -1 , 
A b » 3 • 10- 10 GeV -1 , 

A T « 4 • 1(T 12 GeV" 1 . (17) 

for m susy ~ 100 GeV. There are other loop contributions to neutrino masses and 
a complete discussion is given in ]/], 13, 18|. We have retained here only those 



contributions which are known [|18[ to be dominant in case of mSUGRA and 
GMSB. The additional contributions not included in the text come from, (a) R 
parity violating mixing of the charged sleptons with Higgs fields (b) sneutrino 
exchange diagrams through R-parity violating sneutrino-Higgs mixing and (c) 
loop contribution to the tree level neutrino neutralino mixing. These contri- 
butions are sub-dominant as long as the parameters Amf,ABi are suppressed 



1 18]. Such suppression is required purely from the phenomenological point as 
we argue below. It is then consistent to omit these sub-dominant terms for 
the analytical discussion that follows. We however discuss these additional 
contributions in the appendix B. 

The total neutrino mass matrix is given by 

(Mtatiij = Aoeiejkikj + e l e j (A b + A T (1 - 6 a ) (1 - 6 j3 )) . (18) 

The desired hierarchy among neutrino masses is automatically built in the above 
equations in view of typical numerical values of the parameters Aqi t . The tree 
contribution dominates over the rest (unless ki are enormously suppressed) but 
it leads to only one massive neutrino. Switching on the b-quark contribution 
gives mass to the other neutrino, one neutrino still remaining massless at this 
stage. The latter obtains its mass from somewhat less dominant contribution 
due to A T . Note that hierarchy among the first two neutrino masses need not 
be very strong due to similar magnitudes of A bT . The above statements are 
made explicit below which also contains discussion on neutrino mixing. 
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2.1 Neutrino masses and mixing 



The tree- level neutrino mass matrix can be easily diagonalized: 

MtreeUo = diag{0,Q,m U3 } , 

where 



I c 2 s 2 c 3 s 2 s 3 ^ 

-s 2 c 2 c 3 C 2 S 3 
V -S3 C 3 / 



with S2,3 = sin 6 2t 3 and 

tan 02 = ei^i/e2&2 ; tan #3 



(19) 



(20) 



(21) 



The total mass matrix eq. (|T^) , assumes the following form in basis with diagonal 
tree mass matrix: 



U Mtot Ul 



where 



with 



/ a\(A b + A T ) ai (A b a 2 + A T b 2 ) ai{A b a 3 + A T b 3 ) \ 

ai(A b a 2 + A T b 2 ) A b a 2 + A T b\ A b a 2 a 3 + A T b 2 b 3 
\ a x (A b a 3 + A T b 3 ) A b a 2 a 3 + A T b 2 b 3 A uj 2 + A b a% + A T b\ J 

(22) 



a 2 

h 
b 3 



ei£2 

C3 



(&! - k 2 ), 



U)±UJ 

1 

£3^3 



{elhih - k 3 ) + e 2 2 k 2 (k 2 -k 3 )), 



— (elh + e^) 



u = {e\kl + elkl + elklYI\ 
U± = {eikl + elklfl 2 . 



(23) 



(24) 



The subsequent diagonalization can be approximately done if we neglect terms 
ofO(^f). Let 



( ci si \ 
— s\ ci 

V l l 



(25) 
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where 



We then have 



tan 26 



2ai {A b a 2 + A T b 2 ) 



A h (al 



if) + A 7 



U^oMtotUSu? 



m 



\ 








in 



;/ 2 



o \ 








+ 0( 



Ab,T- 

A n 



(26) 



(27) 



/ 



The eigenvalues are approximately given by 

af(a2 



M 2 



(a? + «l) 



??? 



''2 



4, (a? + ai), 



A UJ 2 . 



(28) 



The mixing among neutrinos is described by 

/ 

u = u% ul = 



V 



c 2 ci - s 1 s 2 c 3 
-s 2 ci - SiC 2 C 3 
S1S3 



(29) 



C 2 Sl + C1S2C3 S 2 «3 
-S 2 Sl + C1C2C3 C 2 S 3 

-C1S3 c 3 ) 

Let us now discuss consequences of the above algebraic results. 

(1) It follows from eqs.(17,28|) that the neutrino masses obey the desired hier 

archy: 

< 



m 



m 



l> 2 



<C m. 



u 3 ■ 



(2) The neutrino masses relevant for the solar and atmospheric scales are re- 
spectively given by Aoe 2 k 2 and Abe 2 leading to 



A, 



A 



at in 



Ab 
A 



1 



where e, k represent typical values of e^, k%. It follows that the ratio of the 
solar to atmospheric scales is independent of the R violating parameters and 
depends upon the values of the soft parameters represented by k. One typically 
needs 



KT 1 GeV ; k ~ 10" 



10" 



(30) 



in order to reproduce the scales correctly. This shows in particular that irre- 
spective of details of the SUSY breaking the Higgs-slepton universality (corre- 
sponding to very small values of k ) is unavoidable in this model if neutrino 
masses are to be correctly reproduced. 
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(3) A particularly interesting limit would be the case where k ~ 0(1). This is 
what one expects in a truly non-universal regime of soft masses |19|] . From the 
above, we see that in this limit, one can infact choose e, ~ O(MeV) and generate 
neutrino masses of 0(eV). However only one neutrino would be massive in this 
scenario, with the the other neutrinos being with extremely negligible masses 
and without any relevance for phenomenology. 

(4) If exact flavour universality were to hold between the first two generations 
then ki = k 2 (see eq.(|ll])). In this case a\ as defined in eq.(p3|) would be zero 
leading to s\ = in eq.(^) . The s± is required to be large in order to obtain the 
large mixing angle solution and obtaining this solution would need very sizable 
departures from the flavour universality among the first two generations. We 
quantify these remarks in the next section. 

3 Neutrino mixing and departure from flavour universality 

We derived approximate expressions for the neutrino masses and mixing with- 
out any specific assumption on the soft symmetry breaking sector. The entire 
neutrino spectrum can be parameterized in terms of three and three ki of 
which ki depend upon the soft SUSY breaking parameters. We now quantify 
the amount of flavour universality violations needed for obtaining the most 
preferred large angle solution to the solar neutrino problem. The following two 
parameters are introduced as a measure of universality violation: 

x = (A* - k 3 )/(h + k 3 ) ; y = (k 1 -k 2 )/(k 1 + k 2 ) (31) 

We regard x and y as independent parameters but restrict their variation to 
values between (-1,1) in the numerical analysis that follows. 

The neutrino mixing is determined by the matrix U in eq. ( ^9[ ) . Due to hier- 
archical mass spectrum, the survival probabilities for the solar and atmospheric 
neutrinos approximately assume two generation form. The corresponding mix- 
ing angles 9q and # a t m are given in terms of elements of the mixing matrix U 
as follows: 

sm 2 26 atm ~ 4 UUl - U^f « 0.8 - 1.0 
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sin 2 26> ~ 4 U^U^ ~ 0.75 - 1.0 
sin 2 Ochooz ~ ^e 3 < 0.01 , (32) 

where numbers on the RHS correspond to the required values for these param- 
eters based on two generation analysis of the experimental data H. 

We can convert the above restrictions on #0,0 a tm to restrictions on the 
mixing angles si,2,3 entering the definition of U. The CHOOZ result requires 
l s 2 s 3| < 0.1 and the nearly maximal atmospheric mixing is obtained with 
IC2S3I « -j=. This requires small s 2 and large S3. The solar neutrino mix- 
ing angle defined in eq.(|32]) coincides with s\ in this limit. We thus need 
sin 2 ~ 0.75 — 1. Large value of s\ in turn needs sizable departure from 
flavour universality as will be argued in the last subsection. 

The expressions for mixing angles and masses obtained in the last section 
can be used to approximately determine the allowed ranges of parameters fcj, 
which explain the solar and atmospheric neutrino anomalies. We approximately 
need | S2 1 < V%U e 3 and \s^\ ~ This implies: 

e 2^2 ~ e 3^3 ) 

|eifci| » y/2\U e3 e 2 k 2 \. (33) 

The magnitude of €3/^3 is then approximately fixed by the atmospheric mass 
scale: 

elkl « ^ « ^ , (34) 



2A n 2A 





while the solar scale and mixing angle determines e 2 : 

2 „ ( \ (l + *) 2 (l-y) 2 

Eqs. (|33| , ^4l , ^5| ) allow us to express magnitudes of all ej, fcj in terms of x,y , 
approximately known ^4q,6 an d the experimentally measurable quantities. 

The solar neutrino mixing angle following from eq.(|26|) is given in the limit 
A T < A b by 

tan 2 6 Q « tan 2 ft « ~ ^ . ( 36 ) 

(x - yY 



We have used eq.(33) in deriving the above relation. It is clear that large Q\ 



require sizable departure from flavour universality, i.e. sizable y. Moreover, 
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one typically needs \x — y\ « 2|J7 e 3y(l — x)| in order to obtain a sizable solar 
neutrino mixing angle. 

We now numerically determine the region in the x, y plane needed to re- 
produce the required ranges in mixing angle and masses. We make use of 
eqs.( |33| - |35|) to determine the approximate input values of 6j, ki in terms of the 
A Q ,A aim and x,y. We allow input values to vary by varying A ,A atm over 
the experimentally allowed ranges. We also randomly choose x,y between -1 
and 1. Through this procedure, we choose a set of 1.5 x 10 5 different values for 
the input parameters e{ , ki . Then we numerically diagonalize the total neutrino 
mass matrix, eq.(|i1f) for each of these values of ej, ki and determine a set of x, y 
values which correctly reproduces the allowed ranges of the solar and atmo- 
spheric neutrino parameters and lead to \U e ^\ < 0.1. We obtain about 2024 x, y 
values leading to the correct description of neutrino anomalies. These points 
in the x,y plane are displayed in Fig.(l). This figure, based on complete diag- 
onalization clearly shows the features obtained through approximate formulas. 
All the allowed values of x and y are in the range —0.9 0.6 and sizable de- 
parture from universality is clearly seen. Also most points satisfy approximate 
equality \x — y\ ~ 2U e sy needed to obtain large solar neutrino mixing angle. As 
an illustration, we give below a typical set of e«, ki which correctly reproduces 
all the parameters: 



e 3 ~ 0.1 GeV 
e 2 ~ 0.031 GeV 
ei ~ 0.087 GeV 



k 3 ~ 1.1 • 10~ 3 

k 2 ~ 3.5 • 10~ 3 (37) 
fci ~ 9.1 • 10~ 4 



Typically, one needs ej ~ O(10 _1 GeV) and ki ~ 10~ 3 as argued before. 

Let us now compare above phenomenological restrictions with expectations 
based on specific framework like mSUGRA. In order to obtain correct neutrino 



masses one needs parameters ki (11) to be suppressed, typically k ~ 10~ 3 — 
10~ 4 as in eq. (|30|) . The other constraint is that y should be 0(1). The k{ 
provide a measure of the Higgs-slepton universality violation. Typical value of 
ki obtained in mSUGRA follows from eq.(l) and is in the range required from 
phenomenology. Thus mSUGRA provides a very good framework to understand 
neutrino mass hierarchy as has been demonstrated in number of papers through 
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detailed numerical calculations J?], |n], [l3| . However mSUGRA would not be able 
to provide the required value of y. This can be seen as follows. Theoretically, 
y can be approximately written using eq.(|TT|) as follows: 



utan^(B 1 - B 2 ) - (m 2 ~ - m 2 -, ) 

v ~ _ 1 v ±k (<va\ 

y /xtan/3(A5i + AB 2 ) - (Am 2 + Am 2 ,) ' V ' 

where we have neglected terms of order (Am?) 2 , (ABi) 2 etc. Within mSUGRA, 
y is identically zero at the high scale as B\ = B 2 and m 2 , 1 = m 2 2 due to the 
universal boundary conditions. At the weak scale, this universality condition is 
broken solely by RG evolution. In the limit of neglecting first two generation 
Yukawa couplings, y is identically zero even at the weak scale. A rough estimate 
of parameters appearing in y can be obtained by approximately integrating the 
RG equations, eqs.(42,|43|) given in appendix A. We see that 



m~, — m 2 -, i / m 

1 2 ~ - | — - I « 10" 



Am 2 + Am 2 6 \ 

V ^ - 10~ 4 . (39) 



(gl ~ B 2 ) ^Ifm,^ 2 
ABi + AB 2 ~ 6 V m b . 



Together they would imply very small value for y ~ instead of the required 
value of 0(1). Thus universal boundary conditions of mSUGRA cannot lead to 
a large mixing angle solution to the solar neutrino problem. 

One need not consider a generic non-universal scenario to introduce univer- 
sality violations. It is clear from the forgoing discussion that one only needs 
small Higgs-slepton universality violation as well as flavour violation of simi- 
lar magnitude to obtain a large solar neutrino mixing angle. These violations 
can come from either non-universal slepton mass terms or from non-universal 
B-terms or both. It is possible that such a scenario can arise from a higher 
theory of flavour either based on string theory [20| or through abelian [21 



or 



non-abelian [22] flavour symmetries. However, for a phenomenological under- 
standing it is clear that the existence of Higgs-slepton and flavour universality 
violations at the high scale would lead to correct neutrino mass spectrum at the 
weak scale. Knowing the value of x and y required for a correct neutrino spec- 
trum at the weak scale, it is possible to estimate the amount of non-universality 
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required at the high scale. For example, using y, we have in the limit of ne- 
glecting contributions from AB terms, the required slepton flavour universality 
violations to be of order: 

to 2 -, (0) m 2 ~, (0) 

m \ (0) - ml, (0) , g K, g (0) + ^(0)) + 2 y - ^ (40) 

where ^m 2 ^ represents the correction to the high scale Higgs mass due to RG 
scaling. From the above we see that for a large negative y, to-^O) should be at 
least a factor of 3 times larger than to| 2 (0). Introducing such a non-universality 
at the high scale would lead to the correct neutrino spectrum at the weak scale. 

A similar analysis can also be considered for the B-terms ||. Moreover, 
such a pattern can be incorporated naturally in non-minimal models of GMSB 
1 23 ] . In these models, identical gauge quantum numbers of all sneutrinos assure 



almost universal sneutrino masses at the weak scale as in the case of mSUGRA. 
In contrast, there is no natural reason within these models for the flavour uni- 
versal B parameters. In fact, the B parameters are assumed to vanish in the 



minimal version of the scheme [24, 25 1. Thus the universality of B parame- 



ters at supersymmetry breaking scale holds by default. It is possible to choose 
non-universal and non-zero B\p terms to start with in this model. This does 
not significantly influence the conventional phenomenology of the minimal ver- 
sion as long as the parameters are much smaller than the /i-parameter in 
the superpotential. But it allows the LMA solution as has been demonstrated 
through a detailed numerical work ||. 

From the above discussion we see that the phenomenological requirement of 
k ~ 0(1O~ 3 ) and x,y ~ O(l) leads to a specific class of non-universality at the 
high scale. A generic non-universal soft spectrum might lead to much larger 
class of solutions. Such an analysis is beyond the scope of present work. 

4 Comments 

Supersymmetric model with bilinear R parity violations provides a potentially 
interesting framework to study neutrino masses and mixing. The dominant 
sources of neutrino masses can be parameterized in this scenario in terms of 
three dimensionful parameters and three dimensional parameters k{. The ki 
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depend on the structure of soft supersymmetry breaking terms at the weak scale. 
We have tried to obtain phenomenological restrictions on £j and fc, without 
making specific assumptions on the values of the soft supersymmetry breaking 
parameters. While neutrino masses can be suppressed by lowering the overall 
scale 6j of R parity violation, phenomenologically preferred hierarchy in neutrino 
masses require that both ej and fc, are suppressed, see eq. (|30|) . fc, provide 
a measure of the Higgs-slepton universality and suppression in their values 
indicate very small amount of this violation. Such violation of universality is 
already built in the mSUGRA and GMSB scenario. 

A large solar neutrino mixing angle can be obtained consistently within 
these scenarios only if flavour universality violations in the soft parameters of 
the first two generations are almost as large as the violation of Higgs-slepton 
universality. This feature does not emerge in models where these universality 
violations are generated solely by RG scaling as in the case of mSUGRA. Thus 
mSUGRA seems more suitable to describe the less preferred small mixing angle 
solution to the solar neutrino problem. 

We concentrated throughout on the most dominant sources of neutrino 
masses in this theory. This is a good assumption in case of small universal- 
ity violation. The other sources of neutrino masses would become important in 
case of large universality violation. It is not unlikely that these contributions 
could also lead to a large solar neutrino mixing angle in such scenarios. 

Another way to achieve the Large Mixing Angle solution is to consider 
some [^, |2^] or all of the dimensionless lepton number violating couplings 
(A, A' ) to be present in the superpotential. Several features of the neutrino 
mass spectrum like hierarchy and large mixing are still preserved within these 
models making them phenomenologically viable. However, unlike the bilinear 
model considered in this work, these models are less constrained simply due to 
the large number of additional couplings present in the model. This can lead 
to the large solar neutrino mixing angle even without relaxing the universality 



constraints on the soft spectrum 
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5 Appendix A 



The renormalization group equations for various parameters appearing in the 
soft scalar potential are basis dependent. We have chosen a specific basis in 
which bilinear terms in the potential are kept in the superpotential till the weak 
scale. These terms are rotated only after evolving to weak scale. We collect 
here RG equations for relevant parameters with this specific choice. They differ 
for example from the ones derived in [18| where relevant rotation is performed 



at each scale. The following equations follow in a straightforward manner from 



the formalism given by Falck[29]: 



J t <>H< = rni, H ,(--Y»--Y T --Y b ), (41) 

^N) = 3Y b (m 2 Q3 + m 2 D3 +m 2 H ,+A 2 b ) 

- Yf(ml+m 2 Ei +m 2 H ,+Af 2 ), (42) 

j t (^Bi) = A T Y T + 3Y b A b -3YfAf. (43) 

In the above, we have used standard notation for all the soft parameters appear- 
ing in the equations with the exception of A which represents the soft trilinear 
couplings. 

6 Appendix B 

In this appendix we justify the neglect of additional contributions to neutrino 
masses not included in the main text. We also discuss flavour violating processes 
/x —* ej and show that the corresponding branching ratio is very small in the 
present context. 

Detailed analysis of the additional 1-loop diagrams contributing to neutrino 
mass matrix has been done in 13, |18|| . While Refs. |7[ calculate all 



the 1-loop self-energy diagrams to the 7x7 neutrino-neutralino mass matrix 
and re-diagonalise it, Ref.JjO]] follows the effective mixing matrix approach. In 
addition to the contributions considered in the text, large contributions are also 
expected from diagrams which are not Yukawa suppressed, thus involving only 
gauge vertices. These can be visualized as diagrams with two R-parity violating 
mass insertions proportional to Am?, ABi as given in eq.(^|), with neutralino 
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(chargino), sneutrino (charged slepton) and neutral Higgs (charged Higgs) in 
the loops (l^, [L8|] . Typical magnitude of these diagrams is given by 

g 2 

Mij » -^2^jK k 'j m ^usy ( 44 ) 



with 

fa* 



ciAmf + c 2 ABf 

9 

mf~. 



m susy is a typical supersymmetry breaking scale and ci 5 2 are coefficients of 
order one following from the scalar mass matrices of the model. k\ are similar 



to parameters k{ defined in eq.(ll). It is natural then to choose k\ ~ -j^ki for 
order of magnitude estimates. Comparing the 1-loop gaugino contribution with 
the 6-quark contribution (eq.(|l^)) Ai b we obtain 

M v„ a 2 /M 2 hkj 

M% l§iPA h \vj m susy V ; 



The numerical value of A b is given in eq.(15). As argued above, we typically 
need h ~ 10 -3 = 1CT 4 . It is seen that the b — quark contribution retained in 
the main text dominates over the gaugino contribution in this case and it is 
consistent to neglect the latter. The other contributions to neutrino masses 
are even less dominant than the gaugino contribution^]. They come from 1- 
loop diagrams with two Yukawa vertices. These can be seen as a) diagrams 
with A and h T vertices with a R-parity violating mass insertion in the internal 
line connecting charged slepton and charged Higgs, and b) diagrams with h T 
couplings at both the vertices with two R-parity violating mass insertions pro- 
portional to the sneutrino vev. Both these sets of diagrams are suppressed by 
the r- Yukawa coupling. They have been analyzed in detail in Ref.^] where it 
has been shown that they can become comparable in magnitude to A T in large 
tan j3 regions. However as we have seen earlier this contribution is always sub- 
dominant compared to the contribution from bottom Yukawa couplings, Aj,. 
Thus it is justified to neglect these contributions within the present analysis. 
Effects of Basis Rotation up to higher order in e : We now generalize the 

basis (|5|) to higher order in e and discuss its consequences. Such generalization 
3 For a detailed discussion of the various diagrams in mass insertion approximation, see 
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becomes necessary for discussion of flavour violating transitions such as /i — > cy. 
Eq.(||) can be re-rewritten as follows: 



I H X \ ( 

Li 
L 2 

V L 3 / 



1 -2 



V 



1-3* 

-ei 



_i;2 _ 

2 £ 1 

\ hh 1 

■5 hh 



(2 

1 -2 



£3 



2 e l £ 3 
1 :2 



4 hh i - 1 ^ y 



V ^3 / 



,+0(6 3 ) 

(46) 



where = (e^)/ '/i, e 2 = + e| + e§. The V^o/j in eq.(f|) assumes the form 

V soft = m 2 Hi \H$\ 2 + m 2 H2 \H° 2 \ 2 + m%\n\ 2 + Am 2 h (vi*H° + c.c) + (-fiBH^H^ + 
- ejASj (viH® + c.c) - - ^2 htji.Arri 2 + Am 2 )(vi*i)j + c.c.) , 



i<7 



where 



2 2 /-i ~2\ , 2 -2 



m, 7 . = m 



2 ~2 , 2/i ~2\ 
H' e i + m uA 1 ~ e i)' 



B = B fl (l-e 2 )+B i e 2 . 



(48) 



The rotation has generated off-diagonal flavour violating sneutrino mixing terms 
at 0(e 2 ). Since these terms conserve lepton number, they do not directly con- 
tribute to the neutrino masses but lead to flavour violating transitions such as 
fx -> ery. 

The rotation in eq.([46|) induces mixing among the charged leptons which 
were diagonal to start with. Define the charged lepton mass matrix as 



then 



where 



U Mi e c , 

1 h\di ei£ 2 h 2 f 2 hhhzh ^ 

Mi = txhhxh h 2 d 2 hhhzh 

\ hhhifi e 2 i 3 h 2 f 2 h 3 d 3 J 



di = vi(l + -{e 2 - |e| 2 ) + e 2 ki - efki , 



(49) 



fi = -vi + h 
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The k{ appearing in above are denned in eq.(pA|) and they signify sneutrino vev 
contribution to the charged lepton mass matrix. As argued in the text, fcj are 
required to be small ~ (10~ 3 — 10 -4 ) in order to account for the correct neutrino 
masses. It then follows that sneutrino vev contribution to each element in Mi 
is suppressed compared to the corresponding contribution of v\. Thus this 
contribution can be neglected while diagonalizing Mi in any realistic theory. 
Even after neglecting it, the C(e 2 ) contribution does produce additional mixing 
among charged leptons that is not Yukawa suppressed. This is easily seen in 
the simplified case of two generation. The 2x2 version of the charged lepton 
mass matrix is obtained from eq.(^) by setting £3 = 0. The following rotation 
on the basis (ei,e2) is needed to diagonalize the charged letpon masses: 




1 \hh 




(50) 



(e, \i) here refers to the flavour basis. This additional rotation affects the neu- 
trino mixing terms in eq.(|47|) which can be re- written in the flavour basis as 



T/ 2 1 ttQ i 2 1 2 1 ttO i 2 1 2i~i2i 2|~|2 

V so ft = mHx\ H l\ + ^hJ-^I + ™4 1 We\ + m i 2 \ u fi\ 

+ (Ampiz^iif + Arr4e 2 r>*Hf + c.c 

- (jiBHiH® + c.c) - (exABxVeltf + e 2 AB 2 u^ + c.c 
1 

: 2 



eie2-{Am\ + Am 2 2 - ml 2 + mlji^D^ + c.c). (51) 



One sees that there are no additional lepton number violating mass terms 
other than present at 0(e). Thus discussion on additional contribution to neu- 
trino masses just given remains unchanged. However, eq.([5l|) contains lepton 
conserving but flavour violating contribution proportional to Pe^V This can 
lead to process such as /j, — > erf. The branching ratio for this process is given 
by 

127T 2 

BR (^ e 7 ) = -^IBI 2 . (52) 

^ F ii 

In the present case, the amplitude B arises due to insertion of the flavour violat- 
ing sneutrino mass term given in the last term in eq. ( |5l| ) . This is approximately 



given by [31] 



151 ~ 16^2^ fc ' (53) 
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where k is a typical magnitude of ki and m~ is sneutrino (mass) 2 . As already 
argued, we need t{ ~ 0.1 GeV and k ~ 10~ 3 . Given this, last equation is seen to 



give very small contribution to BR{[i —* cy) ~ O(10 



-13 



i \k\ 2 ) which makes 



it unobservable in both present [32] and future |33(] experiments. 
Acknowledgements: We gratefully acknowledge useful discussions with E. J. 
Chun which helped in clarifying several issues connected with this work. 
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Figure 1: Allowed values of x and y for which all the neutrino oscillation con- 
straints are satisfied. The input values of parameters are chosen in a way 
described in the text. 
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